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Executive  Summary 

This  document  summarizes  the  work  performed  at 
the  Naval  Civil  Engineering  Laboratory  (NCEL)  and  Solar 
Energy  Research  Institute  (SERI)  under  the  Independent 
Exploratory  Development  (IED)  project“Water  Production 
Dehydrator.”  The  purpose  of  this  effort  was  to  assess  the 
feasibility  of  developing  a  desiccant  system  to  produce 
potable  water  from  atmospheric  sources  that  is  compatible 
with  military  constraints.  The  goals  of  this  effort  were:  (1 ) 
to  examine  desiccant  technology,  investigate  methods  of 
using  available  desiccants  to  collect  atmospheric  moisture, 
(2)  develop  a  conceptual  model  of  a  desiccant  water  produc¬ 
tion  system,  and  (3)  develop  a  mathematical  model  to 
simulate  the  operation  of  the  conceptual  model. 

The  results  show  that  a  desiccant  system  can  pro¬ 
duce  large  quantities  of  potable  water  using  relatively  small 
amounts  of  fuel  for  heat  and  fan  power.  A  solid  desiccant 
based  system,  commercially  available  from  Mitsubishi 
Limited,  has  produced  as  much  as  3.5  gallons  water  per 
gallon  fuel  burned.  Because  it  appeared  that  a  liquid  based 
system  may  have  greater  fuel  efficiency  and  more  compact 
size,  a  hypothetical  system  was  analyzed.  A  comparison  of 
the  maximum  theoretical  performance  of  a  liquid  desiccant 


based  system  versus  a  solid  system  showed  the  liquid 
system  to  have  155  percent  of  the  performance  potential  of 
the  solid  system.  Based  on  this  analysis,  a  conceptual  design 
for  a  liquid  desiccant  system  was  developed.  The  concep¬ 
tual  system  includes  working  fluid  (liquid  desiccant),  ab¬ 
sorber,  heat  exchanger,  distillation  column,  condenser,  and 
subcooler  arranged  as  shown  in  Figure  1 .  This  system  was 
analyzed  using  a  FORTRAN  code  developed  specifically 
for  this  project.  The  code  actually  sized  the  components 
based  on  operating  conditions  and  performance  require¬ 
ments.  The  results  are  shown  in  Table  1  and  are  compared 
to  the  performance  of  the  commercial  unit. 


Table  1 .  Estimated  Performance 


Commercial  Unit  Liquid  Unit 


Water/Fuel  Ratio  3.5  gal 

Power  13.3  kW 

Weight  10,0001b 

Cost  $100,000 


7.3  gal 
6  kW 
4,175  lb 
$33,300* 


♦Components  only 


atrout 


Figure  1 .  Atmospheric  moisture  collection  system  schematic. 
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INTRODUCTION 

Water  is  perhaps  the  most  essential  requirement  for 
human  survival.  Particularly  in  arid  climates,  water  is 
respected  as  a  life  sustaining  essential.  Even  in  the  most 
arid  climates  some  moisture  always  exists  in  the  atmos¬ 
phere  in  the  form  of  humidity,  even  though  no  surface 
water  may  exist.  This  water  is  not  available  directly,  i.e.,  it 
cannot  be  simply  sponged  out  of  the  air  and  squeezed  into  a 
water  glass.  However,  it  is  possible  to  remove  some  of  the 
water  by  condensing  it  out  of  the  air  through  cooling  and/or 
pressurizing.  These  techniques,  when  used  alone,  are  en¬ 
ergy  intensive  and  are,  therefore,  not  of  particular  interest. 
It  has  been  shown  (Ref  1)  that  using  a  desiccant  material, 
substantial  quantities  of  potable  water  can  be  extracted 
from  air  with  absolute  humidity  as  low  as  3  g  HjO/kg. 
Previous  efforts  using  this  technique  have  focused  on  using 
solid  desiccants  such  as  silica  gel  and  molecular  sieves. 

Liquid  desiccants  are  also  effective  at  absorbing 
moisture  and  could  be  used  in  producing  water.  While  the 
basic  process  would  be  the  same  with  solid  or  liquid  desic¬ 
cant,  the  properties  of  a  liquid  desiccant  allow  using  equip¬ 
ment  that  may  result  in  a  more  compact  and  energy  effi¬ 
cient  system  then  possible  when  using  a  solid  desiccant. 

The  focus  of  this  project  was  using  a  liquid  desiccant 
(such  as  triethylene  glycol)  in  an  absorption-distillation 
cycle  much  like  that  used  in  the  petroleum  refining  process. 
The  concept  investigated,  referred  to  as  the  atmospheric 
moisture  collection  system  (AMOCS),  was  developed  into 
a  conceptual  design  and  analyzed  to  predict  its  perform¬ 
ance  potential.  The  performance  predictions  were  based  on 
empirical  data  available  in  literature  (Ref  2,  3,  4,  and  5). 

BACKGROUND 

Humidity  is  water  vapor  that  exists  in  air  as  a  result  of 
the  stoicbastic  nature  of  molecules.  At  a  given  tempera¬ 
ture,  a  portion  of  a  set  of  water  molecules  existing  as  water 
will  have  enough  energy  to  be  liberated  from  the  the  rest  of 
the  liquid  and  take  on  the  gaseous  vapor  form.  Concur¬ 
rently,  a  certain  number  of  water  molecules  existing  at  the 
vapor  state  will  strike  the  liquid  surface  and  give  up  energy 
to  assume  the  liquid  form.  Equilibrium  is  defined  as  when 
the  rate  of  liquid  to  vapor  transition  is  equal  to  the  rate  of 


vapor  to  liquid  transition.  Depending  on  the  temperature 
and  pressure,  the  equilibrium  concentration  of  water  vapor 
above  a  reservoir  of  liquid  water  can  be  calculated.  Thus 
we  have: 

w  =  w(T,  P)  (1) 

where:  w  =  absolute  humidity  (kg  FL,0/kg  air) 

T  =  absolute  temperature  (°K) 

P  =  absolute  pressure  (N/m:) 

The  exact  relationship  between  w,  P,  and  T  is  ex¬ 
tremely  nonlinear,  but  can  be  approximated  or  looked  up  in 
psychrometric  charts  such  as  shown  in  Figure  2. 

Some  of  this  moisture  can  be  condensed  (forced  to 
assume  the  liquid  form)  by  lowering  the  temperature  of  the 
vapor/air  mixture.  Unless  a  low  temperature  heat  sink  is 
available,  however,  this  process  requires  the  expenditure  of 
mechanical  work.  Furthermore,  the  heat  transfer  coefficients 
for  condensation  of  water  vapor  in  the  presence  of  air  arc 
undesirably  low.  If  water  is  to  be  removed  from  air  with 
minimum  work  requirements  and  maximum  heat  transfer 
coefficients,  another  technique  must  be  used. 

One  technique  that  appears  promising  is  a  combina¬ 
tion  of  a  desiccant  debumidifier,  a  desiccant  regenerator, 
and  a  condenser.  A  desiccant  is  a  material  that  has  an 
extremely  high  affinity  for  water  vapor.  The  desiccant  will 
absorb  water  vapor  from  the  air  until  the  desiccant  is 
saturated  with  water.  Desiccant  saturation  occurs  when  the 
desiccant  moisture  density  is  in  equilibrium  with  the  abso¬ 
lute  humidity.  The  saturation  point,  measured  in  mass  of 
water  per  mass  of  dry  desiccant,  is  a  function  of  desiccant 
chemistry,  temperature,  and  pressure.  Thus  wc  have: 

v  =  v  (T,  P)  (2) 

where:  v  -  equilibrium  desiccant  moisture  density 
(ibm  H,0/lbm  desiccant) 

The  exact  relationship  between  temperature,  pres¬ 
sure,  and  equilibrium  desiccant  moisture  density  is  highly 
nonlinear.  Equilibrium  curves  arc  available  in  the  litera¬ 
ture  for  many  desiccants.  Curves  for  tncthylcnc  glycol,  the 
desiccant  chosen  for  this  study,  arc  shown  in  Figure  3. 

Although  the  use  of  desiccants  permits  air  dchydra- 


n  without  the  expenditure  of  mechanical  work  (other 
m  that  required  to  move  the  air),  the  water  is  still  not 
triable  in  liquid  form.  It  has,  however,  been  isolated 
m  the  air.  To  remove  the  water  from  the  saturated  des- 
ant,  the  desiccant  must  be  heated.  This  lowers  the 
jilibrium  moisture  density  of  the  desiccant  and  vapor- 
s  ±z  water  into  steam.  This  process  of  stripping  the 
stccant  of  its  water  requires  heat,  but  no  thermodynamic 
irk.  It  may  appear  now  that  the  processes  of  saturating 
:  desiccant  and  then  stripping  it  have  only  resulted  in  a 
um  to  the  original  state,  plus  the  expenditure  of  some 
it.  The  important  result  is  that  the  water  is  now  in  the 
m  of  steam  rather  than  water  vapor  mixed  with  air.  The 
;her  heat  transfer  coefficients  obtained  when  condensing 
iter  vapor  isolated  from  air  translates  into  much  smaller 
odenser  area  requirements.  Further,  the  steam  tempera- 
e  is  well  in  excess  of  ambient,  so  natural  cooling  may  be 
:d  to  condense  the  steam  to  liquid  water. 

ECHNICAL  APPROACH 

This  project  was  divided  into  three  distinct  phases: 
i  literature  search,  (2)  conceptual  design,  and  (3)  concep- 
il  system  performance  analysis.  Descriptive  details  of 
ch  phase  are  included  in  the  following  paragraphs. 

terature  Search 

The  purpose  of  the  literature  search  was  to  establish 
;  current  state  of  air  dehydration  technology  as  well  as 
/iew  any  prior  work  in  atmospheric  water  collection.  It 
is  immediately  clear  that  desiccants  were  by  far  the  most 
pular  choice  for  air  dehydration.  However,  nearly  all  the 
plications  first  reviewed  were  concerned  with  using  solid 
siccants  to  debumidify  airstreams.  References  6,  7,  8, 
i  9  discuss  and  analyze  various  techniques  for  using 
lid  desiccants,  primarily  silica  gel,  fordehumidifying  an 
stream.  Typically,  these  applications  are  part  of  a  solar 
wered  air-conditioning  system.  The  solar  energy  was 
:d  directly  to  heat  and  dry  (regenerate)  the  desiccant.  A 
iple  analysis  indicates  that  the  collector  size  for  solar 
jeneration  of  desiccant  in  a  water  production  system 
iuld  be  unpractical  for  a  mobile  system.  This  is  due  to 
:  high  latent  heat  of  vaporization  of  water  when  con- 
tsed  on  a  desiccant,  the  relatively  low  power  level  of 
ar  insolation,  and  the  need  to  produce  relatively  large 
intities  of  water.  Therefore,  these  references  were  of  pe- 
heral  interest. 

Liquid  desiccant  based  air-conditioning  systems  were 
o  reviewed.  References  10  and  11  discuss  various 
hniques  for  regenerating  liquid  desiccants  such  as  litfa- 
n  chloride  and  triethylene  glycol.  The  major  thrust  of 


these  analyses  was  to  regenerate  the  desiccant  without 
regard  to  collecting  the  moisture  (as  was  that  of  the  solid 
desiccant  research).  These  references  were  also  of  only 
limited  interest. 

The  natural  gas  industry  has  been  using  various  liq¬ 
uid  desiccants  to  dehydrate  natural  gas  for  many  years. 
This  application  was  clearly  the  most  directly  applicable  to 
the  current  project.  A  thorough  review  of  the  industry's 
publications  revealed  a  weal’h  of  knowledge  in  the  form  of 
empirical  desiccant  regeneration  and  absorption  data  (Ref 
2, 3, 4,  and  5)  and  design  data.  Reference  12  also  provided 
industry  standard  design  guidelines  for  direct  contact  gas 
coolers  and  condensers. 

In  the  course  of  the  literature  search,  a  commercial 
device  called  the  Mobile  Water  Collection  System  (MWCS) 
built  by  Mitsubishi  Electric,  Ltd  was  identified.  A  review 
of  MWCS  literature  showed  that  it  was  based  on  a  solid 
desiccant  operating  in  a  batch  mode.  Its  performance  is 
quoted  as  3.5  gallons  water  per  gallon  of  fuel  burned  and 
400  gallons  of  water  produced  in  a  24-hour  period.  Thus, 
the  feasibility  of  producing  water  from  atmospheric  sources 
is  demonstrated.  The  goal,  therefore,  of  the  conceptual 
system  design  for  AMOCS  was  to  determine  whether  there 
was  potential  capability  to  perform  substantially  better 
than  the  commercial  system. 

Conceptual  Design 

The  original  concept  of  the  AMOCS  was  based  on  a 
solar  regenerated  solid  desiccant  As  described  previously, 
solar  regeneration  is  impractical  for  this  type  of  applica¬ 
tion.  Further,  because  the  liquid  desiccant  systems  had 
been  so  successfully  used  in  the  natural  gas  industry,  liquid 
desiccant  was  selected  as  the  basis  for  the  conceptual 
design.  Finally,  the  substantial  performance  improvement 
being  sought  required  a  departure  from  previous  concepts. 
Because  the  commercial  MWCS  was  based  on  a  solid  des¬ 
iccant,  and  there  was  no  evidence  that  liquid  desiccant  sys¬ 
tems  had  been  considered,  a  liquid  based  system  approach 
was  selected. 

An  AMOCS  can  operate  in  either  a  continuous  or 
batch  process  mode.  The  solid  based  MWCS  operates  in  a 
batch  process  mode  only.  While  one  adsorption  chamber  is 
being  regenerated  and  producing  steam  for  condensation, 
the  other  chamber  is  being  saturated.  Once  saturation  and 
regeneration  are  complete,  the  processes  are  switched.  The 
chamber  that  has  been  fully  regenerated  is  now  switched  to 
an  adsorption  mode  and  the  fully  saturated  chamber  is 
regenerated.  The  system  schematic  is  shown  in  Figure  4. 
Similarly,  a  liquid  desiccant  based  system  can  operate  in  a 
batch  process  mode  or  a  continuous  mode.  The  continuous 
pmcess  mode  was  selected  for  the  conceptual  design  based 
on  the  experience  of  the  natural  gas  industry. 


A  continuous  process  liquid  desiccant  based  AMOCS 
requires  five  major  components.  These  components  are 
shown  in  the  schematic  of  the  conceptual  design  in  Figure 
1 .  The  absorber  exposes  the  lean  (water  poor)  desiccant  to 
the  ambient  air.  Some  of  the  moisture  in  the  air  (hu  tiidity) 
is  absorbed  by  the  desiccant  and  its  latent  heat  of  vaporiza¬ 
tion  is  released.  The  desiccant,  now  rich  (water  rich),  is 
pumped  through  the  heat  exchanger  where  it  receives  heat 
from  the  lean  desiccant  on  its  way  to  the  absorber.  From 
the  heat  exchanger,  the  rich  desiccant  is  fed  to  the  regen¬ 
erator  where  it  is  heated  and  stripped  of  most  of  its  water. 
The  now  lean  desiccant  is  then  fed  back  through  the  heat 
exchanger  and  on  to  the  absorber.  This  constitutes  the 
continuous  process  water  collection  cycle.  The  water 
stripped  from  the  desiccant  in  the  regenerator  is  now  in  the 
form  of  steam.  The  steam  is  directed  to  a  condenser  to 
liquify  the  steam  using  forced  air  cooling.  The  condensate 
leaving  the  condenser  is  still  quite  hot  ( =  100°C)  and  is 
cooled  further  by  a  subcooler  also  using  forced  air  cooling. 

Conceptual  System  Performance 
Analysis 

Selecting  the  components  for  the  conceptual  design 
was  based  on  natural  gas  industry  standards.  Several 
choices  were  available  for  the  absorber  and  regenerator. 
Baffle  tray  columns,  spray  chambers,  packed  columns,  and 
pipeline  contactors  are  valid  options.  Because  of  the  vol¬ 
ume  of  air  that  must  be  exposed  to  the  desiccant  in  the 
absorption  process  and  the  tendency  for  sprays  to  become 
entrained,  the  spray  chamber  and  pipeline  contactor  op¬ 
tions  were  eliminated.  Packed  columns  were  selected  over 
baffle  tray  columns  due  to  their  higher  spatial  efficiency. 
A  packed  column  was  selected  for  the  regenerator  for 
similar  reasons.  The  heat  exchanger  was  a  simple  plate 
counter  flow  device.  The  condenser  and  subcooler  were 
simple  air  cooled  heat  exchangers. 

Desiccant  Selecting  the  desiccant  for  the  concep¬ 
tual  design  required  reviewing  several  compounds  cur¬ 
rently  used  in  industry.  Of  particular  concern  was  the 
toxicity  of  the  desiccant  to  humans  upon  ingestion.  It  was 
anticipated  that  trace  quantities  of  the  desiccant  would  be 
left  in  the  condensate  after  distillation.  Tnethylene  glycol 
was  selected  based  on  the  availability  of  hygroscopic  data 
and  low  toxicity  characteristics. 

To  simplify  later  analysis,  a  design  point  was  se¬ 
lected  and  a  conceptual  design  was  developed  around  that 
point.  The  design  conditions  were: 

Tot  -  30°C 

RH  -  30% 


(w  »  0.0078  kg  H-,0/kg  air) 

Output  ■=  10  gallons  per  hour 

where:  TM  -  Temperature  dry  bulb 
RH  =  Relative  humidity 

Distillation  Columns.  Once  the  component  types 
had  been  specified  and  the  design  point  selected,  the  com¬ 
ponents  could  be  sized  and  priced.  The  distillation  column 
requirements  could  be  analyzed  using  either  a  continuous 
exchange  model  or  an  equilibrium  tray  model.  The  equi¬ 
librium  tray  technique  was  selected  due  to  the  simplicity  of 
the  method  and  the  lack  of  requirement  for  mass  and  heat 
transfer  coefficient  data.  A  computer  routine  based  on  the 
equilibrium  tray  approach  using  the  McCabe-Thiele  gra¬ 
phical  method  was  used  to  model  the  distillation  column. 
The  distillation  model  analysis  is  shown  in  Appendix  A. 
The  program  "DIST5.FOR."  written  in  Microsoft  FOR¬ 
TRAN  is  also  listed  in  Appendix  A. 

The  final  distillation  column  design  values  were: 

Theoretical  Trays  »  3 

Equiv.  Height  Packing  -  28  in.  (Flexipack  1 Y) 

Column  Height  =  84  in. 

Column  Diameter  =  6  in. 

Pressure  Drop  =  0.84  in.  H^O 

Cost  =  $4,000  (including 

packing) 

The  other  parameters  associated  with  the  distillation 
column  that  bad  to  be  calculated  were  the  distillation  peak 
temperature,  reflux  ratio,  and  the  reboiler  load.  The  peak 
allowable  temperature  was  375°F.  the  maximum  tempera¬ 
ture  of  triethylene  glycol  before  disassociauon  occurs.  This 
maximum  distillation  temperature  also  restricts  the  level  of 
regeneration  that  can  be  achieved  at  atmospheric  pressure. 
The  maximum  purity  of  desiccant  leaving  the  distillation 
column  was,  therefore,  calculated  to  be  98  percent  (2 
percent  water).  The  purity  of  the  steam  leaving  the  distilla¬ 
tion  is  controlled  by  the  reflux  ratio  (portion  of  condensate 
returned  to  distillation  column  divided  by  the  portion  re¬ 
moved  from  the  desiccant  cycle).  Based  on  the  toxicity 
levels  for  triethylene  glycol,  the  corresponding  water  pu¬ 
rity  level  must  be  greater  than  99.99  percent  (less  than  0.01 
percent  triethylene  glycol).  This  requires  a  reflux  ratio  of 
at  least  0.27. 

Absorption  Column.  The  absorption  column  was 
designed  using  the  same  technique  used  for  the  distillation 
column.  A  separate  computer  routine  was  written  to  model 
the  absorption  column  and  is  listed  in  Appendix  B.  The 
final  design  values  for  the  absorption  column  were: 
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Theoretical  Trays  =  1 

Equiv.  Height  Facking  =  48  in.  (Structured 

Flexipak  Type  2) 

Column  Height  =  48  in.  (external  «  8  feet) 

Column  Diameter  =  48  in.  (external  =  5  feet) 

Pressure  Drop  =  2.0  in.  H 20 

Cost  =  $20,000 

Heat  Exchanger.  The  heat  exchanger  was  de- 
aed  based  on  the  temperatures  of  the  the  saturated,  cool 
iccant  coming  from  the  absorption  column  and  the  hot, 
pped  desiccant  coming  from  the  distillation  column  and 
desiccant  flow  rate.  Assuming  an  overall  heat  transfer 
fficient  of  40  Btu/hr/ ft2,  approximately  1,050  ft2  of 
hanger  area  are  required.  A  plate  type  heat  exchanger 
neet  this  requirement  was: 

Size  =  36  x  24  x  24  in. 

Weight  =  400  lb 
Cost  =  $3,000 

Condenser.  The  condenser  was  designed  based  on 
incoming  steam  flow  rate  and  the  ambient  (design 
nt)  air  temperature.  Its  heat  load  was  calculated  to  be 
1,000  Btu/hr  and  an  air  cooled  heat  exchanger  was: 

Size  -  66  x  39  x  27  in. 

Weight  -  385  lb 
Cost  = $1,470 

Subcooler.  The  subcooler  was  designed  based  on 
temperature  of  the  incoming  condensate  and  a  desired 
put  temperature  of  %°F.  This  output  temperature  was 
:cted  as  the  lowest  practical  temperature  for  forced  air 
ling.  Its  beat  load  was  calculated  to  be  12,000  Btu/hr 
was: 

Size  -  47  x  22  x  22.5  in. 

Weight  =  156  lb 
Cost  = $650 

:  forced  air  cooler  could  be  eliminated  in  favor  of  a 
ling  tank  if  time  is  available  for  natural  convection 
ling. 

Waste  Heat  During  the  conceptual  design,  it  was 
gested  that  waste  beat  from  a  generator  or  other  engine 
ht  be  available  to  beat  the  desiccant  for  regeneration, 
leteimine  the  feasibility  of  this  concept,  the  heat  avail- 
:  from  engine  exhaust  was  compared  to  the  heat  re- 
ements  of  this  system.  The  exhaust  temperature  from  a 
cal  turbo-charged  diesel  engine  is  approximately  750°F. 
lat  exhaust  stream  is  cooled  to  400°F,  the  engine  pro- 
id  roughly  1,000  Btu/hr/hp.  Based  on  a  need  for 


145.500  Btu/hr  in  the  reboiler,  a  150-bp  engint  should 
supply  enough  heat  to  regenerate  the  desiccant.  Additional 
heat  may  be  available  from  the  engine  cooling  jacket. 

RESULTS 

Once  all  the  components  had  been  sized  and  se¬ 
lected,  a  system  analysis  could  be  performed.  Based  on  the 
design  point  operation,  the  system  was  calculated  to  per¬ 
form  as: 

Water  Output  =  10.7  gph 
Power  Requirements 

Desiccant  Loss  =  0.4  gal  1  day  (24  hr) 

Reboiler  =  1 .5  gal  fuel/hr 

Absorption  Column  Fan  =  3.0  kW 

Condenser  Fan  =  1.6  kW 

Subcooler  Fan  =  0.83  kW 

Desiccant  Pump  =0.10kW 

Total  =  5.5  kW  +  1 .5  gal  fuel/hr 

The  conceptual  system  cost  and  weight  breakdown  is: 


Item 

Weight 

(lb) 

Cost 

(S) 

Distillation  Column 

250 

4.000 

Absorption  Column 

1,500 

20.000 

Reboiier 

1.400 

3.670 

Heat  Exchanger 

400 

3.000 

Condenser 

385 

1,470 

Subcooler 

160 

650 

Desiccant  Pumps  (2) 

80* 

510* 

Total 

4,175 

33.300 

•Each  pump  weighs  40  pounds  and  costs  $205.00 

At  this  point  in  the  analysis,  it  is  possible  to  compute 
the  effective  ratio  of  water  produced  to  fuel  consumed 
Assuming  a  heating  value  of  135,000  Btu  per  gallon  of  fuel 
and  20  percent  heat  loss,  the  ratio  is  7.3  gallons  H,0 
produced  per  1  gallon  of  fuel  burned.  Compared  to  the 
ratio  for  the  commercially  available  MWCS  solid  desic¬ 
cant  system,  this  result  is  extremely  favorable  even  though 
it  is  based  on  relatively  conservative  assumptions.  More¬ 
over,  the  electric  power  requirement  of  this  theoretical 
liquid  desiccant  system  is  less  than  half  that  of  the  equiva¬ 
lent  size  MWCS.  Rough  estimates  of  size  and  weight  of 
the  liquid  system  suggest  that  a  10-gallon  per  hour  device 
could  be  built  within  a  8-  by  8-  by  10- foot  frame  (HALF- 
CON)  and  would  weigh  less  than  5,000  pounds.  The  cost 


of  manufacturing  a  liquid  based  system  depends  heavily  on 
the  quantity  produced.  In  1986,  the  estimated  costs  of 
components  for  the  system,  if  purchased  separately,  were 
$33,300. 

Perhaps  the  single  greatest  advantage  of  using  liquid 
desiccant  is  the  greater  potential  water  production  rate  for  a 
given  system  weight  Dividing  the  estimated  water  pro¬ 
duced  per  hour  by  the  system  weight,  we  find  the  specific 
water  productivity  is  2.59  x  10'3  gallon  of  water  per  hour 
per  pound.  The  MWCS  solid  desiccant  system,  on  the 
other  hand,  has  a  specific  productivity  of  5.55  x  10“*  gal¬ 
lons  of  water  per  hour  per  pound.  The  theoretical  liquid 
system  has  nearly  a  five-fold  advantage  in  specific  produc¬ 
tivity. 

The  unique  capability  of  the  desiccant  based  water 
collection  system,  solid  or  liquid  based,  is  its  ability  to 
produce  water  without  a  feed  water  source.  Although  other 
water  purification  systems  such  as  reverse  osmosis  require 
a  feed  water  source,  a  comparison  of  their  performance 
with  a  desiccant  system  might  be  interesting.  An  existing 
reverse  osmosis  system  is  the  Marine  Corps'  600  gallons 
per  hour  (GPH)  Reverse  Osmosis  Water  Purification  Sys¬ 
tem  (ROWPU).  As  its  name  unplies,  it  produces  600  GPH 
of  purified  water  from  salt  water.  It  can  produce  1200  GPH 
from  fresh  water  sources.  The  device  weighs  approxi¬ 
mately  7300  pounds  and  requires  a  30  kW  external  power 
supply.  It  fits  within  a  HALFCON  frame.  Its  specific 
productivity  (system  weight  divided  by  water  production 
rate)  ranges  from  12.2  to  6.1  pounds  per  gallon  of  water  per 
hour.  Clearly,  the  capability  to  produce  water  without  a 
feed  water  source  carries  with  it  a  severe  weight  penalty. 

CONCLUSIONS 

The  effort  documented  here  has  attempted  to  assess 
the  feasibility  of  developing  a  mobile  desiccant  system  to 
produce  potable  water  for  military  use.  There  is  a  commer¬ 
cial  system  based  on  a  solid  desiccant  that  is  capable  of 
producing  3.5  gallons  of  water  per  gallon  of  fuel  burned. 
This  study  was  focused  on  investigation  of  liquid  desiccant 
based  systems.  A  liquid  desiccant  system  was  conceptual¬ 
ized  and  analyzed  using  fundamental  engineering  analysis 
and  computer  models.  The  results  of  the  liquid  system 
analysis  indicate  that  substantially  higher  efficiency  and 
specific  productivity  can  be  achieved  by  the  liquid  system. 

Further  work  is  necessary  before  exploratory  devel¬ 
opment  is  completed.  First,  more  flexible  models  must  be 
developed.  Improved  models  will  permit  examination  of 
the  system  and  the  impact  of  the  following  variables: 

1.  Operation  in  off-design  conditions 


2.  Simulation  of  different  desiccants 

3.  Vacuum  distillation  of  desiccant 

4.  Higher  desiccant  flowTates 

RECOMMENDATIONS 

If  the  improved  computer  models  support  the  indica¬ 
tions  of  the  study  contained  herein,  a  process  model  of  the 
conceptual  system  should  be  developed  and  tested.  As¬ 
suming  the  process  model  further  supports  the  results  of 
this  study,  NCEL  recommends  that  a d  advanced  devel¬ 
opment  model  of  the  conceptual  liquid  based  system  be 
fabricated  and  tested  to  verify  the  performance  predictions. 
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Figure  4.  Mitsubishi  mobile  water  collection  system  schematic  (batch  process). 


Appendix  A 


DISTILLATION  COLUMN  ANALYSIS 


The  distillation  column  is  modeled  using  the  equilib¬ 
rium  tray  approach  (see  Figures  A-l  and  A-2): 


where: 


d  -  molar  flow  rate  of  distillate 
Xd  =  mole  fraction  of  water  in  distillate 
Q  =  condenser  heat 
r  *  reflux  rate 

b  =  molar  flow  rate  of  bottom  product 
Xb  =  mole  fraction  of  water  in  bottom  product 
f  =  molar  flow  rate  of  feed  product 
Xf  =  mole  fraction  of  water  in  feed  product 
-  reboiler  heat 

V  -  total  molar  gas  flow  rate 

L  ”  total  molar  liquid  flow  rate 

Y  *  mole  fraction  of  steam  in  gas 
X  -  mole  fraction  of  water  in  liquid 
P  -  refers  to  any  plate 


The  following  are  known: 

d  **  5.2  LB-moles/hr.  determined  from  water 
requirement  200  gal/day  x  1  day/ 18  hrs  x 
8.34  liygal  ■  92.7  Ih/hr  or  5.2  LB-moles/hr 


SOLUTION  PROCEDURE 

1.  Y,  -  f(X,)  Equilibrium 

V  =  f(X) 

2.  H  =  f(Y) 

3.  Assume  Xp.l.bp.,.  =  f(Xp.|) 

4.  Solve  for  V  and  L„. ,  from  two  simultaneous 

P  P*  I 

equations: 

Vp  Hp  =  Lp-  i  hP- 1  +  M  +  Qc  ener8y  balance 

V  -  L  ,  -  d  mass  balance 

P  9*  I 

5.  Determine  x  .  from 

P  *  1 

V  Y  =  L  X  +  X.d  water  balance 

P  P  p*  1  P*  1  d 

6.  Substitute  X  „ .  for  value  assumed,  iterate  to 

9*  1 

convergence 


Xd  -  0.99992:  from  purity  requirement  7.  Move  to  next  plate  up  until  Xp< ,  /  Xd 

Xu  -  0.145  t98%);  determined  from  maximum 

allowable  boiling  temperature  for  glycol  -  RESULTS  (see  Figure  A-2) 

375°F 


X,  «  0.318  (95%):  based  on  glycol  flow  of  3  gal 
per  lb  H^O  absorbed 

f  -  75.2  LB-moles/hr,  b  -  20.0  LB-moles/hr 
determined  from  simultaneous  solution  of 
total  mass  balance  and  water  balance: 

b  -  f  -  d 
Xfcb  -  X/-X,d 

Q,  *  (d  +  r)  hf| 

Unknowns:  V  ,  L  .,  Y  ,  X  , ., 

p  p*  i  p’  p  ♦  r 

H.h  „N,Q„ 

p*  p  ♦  1  *  ’  ^RJB 


C  This  program  simulates  a  binary  distillation  column 

PROGRAM  DISTILL 
REAL  L(IOO) 

DIMENSION  V(100).  X(100),  X2(I00),  Y(100) 

DIMENSION  HVAP  ( 1 00).  HLIQ(  1 00).  TDIST(  1 00) 

KDFLAG  -  0 

XD  -  .9992 

FEED  -  25.2 

XF  -  .318 

XB  -  .1445 

DIST  -  FEED*(XF  -  XB)/(XD  -  XB) 

STEPCON  -  .001 
HVAPD  -  20709 
HLIQD  -  3242 
ICOUNT  -  1 


A-1 


BETA  =  .05 
1=  1 

X(l)  =  XB 

CALL  EQU1L(X(1).Y(1)) 

RRMIN1  -  (XD  -  Y(1V(XD  -  X(l)) 

RRMIN  -  RRMIN1/(1  -  RRMTN1) 

PRINT.  “ENTER  REFLUX  RATIO:" 

INPUT,  RR 

REFLUX  =  RR*DIST 

QC  =  (DIST+REFLUX)*(HVAPD-HLIQD) 

X(2)  =  .80 

I  FORMAT  (lX,"PLATE",’T",’,X","r,,”V”.”I,';X”) 
PRINT  50 
X)  CONTINUE 

IF  (ICOUNT.GT  2000)  THEN  PRINT,  “NOT  CON- 
ERGING.  TRY  AGAIN  WITH  A  HIGHER 

EFLUX.” 

GOTO  160 
ELSE 

CONTINUE 
END  IF 

IF  (ICOUNT.GT.500)  BETA  -  .2 
IF  (ICOUNT.GT.  1000)  BETA  -  .4 
CALL  SATVAP(Y(I).HVAP(I),TDIST(I)) 

CALL  SATLIQ(X(I+ 1  ),HLIQ(I+ 1 )) 

L(I+1)  -  (QC  -  (HVAP(I)-HLIQ)*DISTV(HVAP(I)  - 
LIQ(I+1)) 

V(I)  -  L(I+1)  +  DIST 

X2(I+1)  =  (V(I)*Y(I)  -  XD’DISTVL(I+1)) 

110  FORMAT  (IX.  “X2  -  “E8.5) 

PRINT  110,  X2(I+1) 

IF  (X2(I+1).GT.0.99)  STEPCON  -  .00001 
IF  (X2(I+1).LE.0.99)  STEPCON  -  .001 
IF  ( ABS(X2(I+1  )-X(I+ 1  ).GT.STEPCON)  THEN 
X(I+1)  =  X(I-t-l)  +  BETA*(X2(I+1)  -  X(I+1)) 
ICOUNT  -  ICOUNT  +  1 
GOTO  100 
ELSE 

CONTINUE 
END  IF 

IF  (X(2).GT.1.0)  THEN 

PRINT,  “X  GREATER  THAN  1.  USE  A 
HIGHER  REFLUX.” 

GOTO  160 
ELSE 

CONTINUE 
END  IF 
ICOUNT  -  1 
BETA  -  .05 
STEPCON  -  .001 

5  FORMAT  (IX.  13,  F10.5,  F8.5,  F6.1.  F7.2,  F7.2) 
PRINT  125, 1.  Y(I),  X(I+1),  TDIST(I).  V(I),  L(I+1) 
IF  (X(I+1  ).GE.XD.AND.KDFLAG.EQ.O)  THEN 
DPLATES  -  I  -  1  +  (XD-X(iy(X(I+l)-X(l)) 


KDFLAG  =  1 
GOTO  130 
ELSE 

CONTINUE 
END  IF 

130  IF  (XD.GT.X(Ifl))  THEN 
I  “  1+1 

CALL  EQUIL(X(I),Y(I)) 

X(I+1)«X(I) 

GOTO  100 
ELSE 

CONTINUE 
END  IF 
PRINT 

150  FORMAT  (IX.  “NUMBER  OF  THEORETICAL 
PLATES  -  ",  F6.2) 

PRINT  150,  DPLATES 
PRINT 
160  END 

SUBROUTINE  EQUIL  (V,S) 

REAL  X(25),Y(25),BY(25),CY(25),DY(25) 

REAL  S,  U.  V,  SEVAL 

CHARACTER  FLAG *4 

INTEGER  I,  N 

DATA  X/0.000, 

.05,.10..15,.20,.25,.30,.35..40..45..50,.55,.60,.65. 

*  .70..75..80..85..875..90..93  3  ,.95  ..975  ..99, 1 .0/ 
DATA  Y/0.0..00024..00072..001 67..00286, 
.00406, .00525. ,00668,  *.00835,.0l026..0l241 , 

.01 456,.01694..01957,.0227..029..044..079,  M05. 

.1 49,.255,.352,.585,.805,1  .(V 

N-25 

IF  (FLAG  .EQ.  ‘DONE-)  GO  TO  2 
CALL  SPLINE  (N,  X,  Y,  BY,  CY.  DY) 

FLAG  -  ‘DONE’ 

2  U-  1.0- V 

S  -  1.0  -  SEVAL  (N,  U,  X,  Y,  BY,  CY.  DY) 
RETURN 
END 
C 

SUBROUTINE  SPLINE  (N,  X,  Y,  BY,  CY.  DY) 

INTEGER  NNM1.IB.I 

REAL  X(N),Y(N),BY(N),CY(N),DY(N) 

REALTY 

NM1-N-1 

IF  (N  .LT.  2)  RETURN 
IF  (N  .TL.  3)  GOTO  50 
DY(1)  =  X(2)-X(l) 

CY(2)  -  (Y(2)-Y(I)yDY(l) 

EX)  101-2,  NM1 
DY(I)  -  X(I+1)-X(I) 

B Y(I)  -  2,*(DY(I-1  )+DY(D) 

CYG+l) «  (Y(I-t-l)-Y (I)VDY ( I ) 

CY(I)  -  CY(I+1)-CY(I) 
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10  CONTINUE 

BY(1)  -  -DY(1) 

BY(N)  —  -DY(N-l) 

CY(1)  =  0. 

CY(N)  -  0. 

IF  (N  .EQ.  3)  GO  TO  15 
CY(1)  =  CY(3V(X(4)-X(2))-CY(2V)X(3)-X(1)) 
CY(N)  -  CY(N-lVX(N)-X(N-2))-CY(N-2V(X(N- 
l)-X(N-3)) 

CY(l) «  CY(1)*DY(1)**2/(X(4)-X(1)) 

CY(N)  -  -CY  (N)*  D  Y(N- 1  )**2/(X(N)-X(N-3 )) 

15  DO  20  I-2.N 

TY  *  DY(I-1)/BY(I-1) 

BY(I)  =  BY(I)-TY*DY(I-1) 

CY(I)  =  CY(I)-TY*CY(I-1) 

20  CONTINUE 

CY(N)  -  CY(N)/BY(N) 

DO  30  IB  -  1.NM1 
I  *  N-IB 

CY(I)  -  (CY(I)-DY(I)*CY(I+1VBY(I) 

30  CONTINUE 

BY(N)  -  (Y(N)-Y(NMl)yDY(NMl)  +  DY(NMl) 
*(CY(NM1)+2.*CY(N)) 

DO  40  I  -  1J4M1 

BY(I)  -  ( Y(I+ 1 )- Y (1)1/ D Y(I)  -  DY(I) 

*(CY  (1+ 1  )+2.*CY(I)) 

DY(I)  -  (CY (1+ 1  )-CY (I)yDY (I) 

CY(I)  -  3.*CY(I) 

40  CONTINUE 

CY(N)  -  3.*CY(N) 

DY(N)»DY(N-l) 

RETURN 

50  BY(1)  -  ( Y (2)- Y ( 1  )y(X(2)-X(  1 )) 

CY(1)«0. 

DY(1)  -  0. 

BY(2)  -  BY(1) 

CY(2)  =  0. 

DY(2)  -  0. 

RETURN 

END 

C 

C 

REAL  FUNCTION  SEVAL  (N,U.X,Y,B,C,D) 
INTEGER  N 

REAL  U.  X(N).  Y(N),  B(N),  C(N).  D(N) 

INTEGER  I,  J,  K 
REAL  DX 
DATA  1/1/ 

IF  (I  .GE.  N)  I  -  1 
EF(U.LT.  X(I))  GO  TO  10 
IF (U  .LE.  X(I+I))  GO  TO  30 
101-  1 
J-N+l 
20  K  -  (I+jy2 


IF  (U  .LT.  X(K))  J  -  K 
IF  (U  .GE.  X(K))  I  =  K 
IF  (J  .GT.  1+1)  GO  TO  20 
30  DX  -  U-X(I) 

SEVAL  -  Y(I)+DX*(B(I)+DX*(C(I)+DX*D(I))) 

RETURN 

END 

SUBROUTINE  SATLIQ  (V.HLIQ.S) 

REAL  X(25),TEMP(25),BT(25),CT(25),DT(25) 
REAL  S,U,V,KEL,HLIQ,A,B,C,D,E,HH20, 
HTEG.HTSOLN 
CHARACTER  FLAG*4 

DATA  X/0.000,.05,.10,.15,.20,.25,.30,.35,.40,.45, 
.50,.55,.60,.65,  *.70..75,.80,.85..875,  .90..933, 
.95,.975,.99,1.0/ 

DATA  TEMP/0.0,.00309,.00926,.02160,.03704, 
.05247..06790..08642, 

*.10802, .13272,.  16049.  .18827..21914..25309. 
.29321,. 34568..41975, 

*.5 1 852..57099..62963..7253 1  ,.78395, .87963 . 

.95062,1.0/ 

N=25 

IF  (FLAG  .EQ.  DONE’)  GO  TO  2 
CALL  SPLINE  (N,  X.  TEMP,  BT,  CT.  DT) 

FLAG  -  ‘DONE’ 

2  U  -  1.0- V 

S  -  324.0  *  SEVAL  (N,  U,  X,  TEMP,  BT.  CT,  DT) 

+  212.0 

KEL  •  (S-32.)*5/9  +  273.15 
DATA  A/5.3 1 409/.B/-.0072001 7/.C/.00001 14153/ 
HH20-  V*(A*(KEL-273.15)  +  B/2.*(KEL**2- 
273.15**2) 

*  +C/3.*(KEL**3-273.15**3))*7.7546 

DATA  D/1.5380E5/.E/5.87E2/ 

HTEG  -  U*(D*(KEL-273.15)  +  E/2*(KEL**2- 
273.15**2))*4.3033E-4 
DATA  HTSOLN/1 549.72/ 

IF  (V  .GT.  0.7)  HTSOLN-O.O 
HLIQ  -  HH20  +  HTEG  -  V*HTSOLN 
RETURN 
END 
C 

SUBROUTINE  SATVAP  (V.HVAP.S) 

REAL  Y(25),TEMP(25),BT(25),CT(25),DT(25) 
REAL  S,U,V,HVAP,KELvA.B,C,D.E.F.G,H.U 

*  H20HTVAPN,TEGHTVAPN,HH20,HTEG.TCRIT 
CHARACTER  FLAG *4 

DATA  Y/0.0,.00024.,00072,.00l67,.00286..00406., 
00525  ,.00668  ..00835..01 026..01 241  ..01 456.. 
01694,.01957,.0227,.029,.044,.079. 
*.105,.I49..255,.352..585,.805,I.(V 

DATA  TEMP/0.0,.00309,.00926,.02160,.03704, 
.05247..06790..08642. 

*.  10802,.  13272,.  16049,.  18827..21914..25309..29321. 


.34568, .41975 

♦.51852..57099..62963..72531,. 78395, 
,87963,.95062,1.0/ 

N=25 

IF  (FLAG  .EQ.  ‘DONE’)  GO  TO  2 
CALL  SPLINE  (N,  Y,  TEMP,  BT,  CT,  DT) 

FLAG  =  ‘DONE’ 

U=  1.0-  V 
FORMAT  (FI  2.3) 

S  =  324.0  *  SEVAL  (N,  U,  Y,  TEMP,  BT,  CT,  DT) 

+  212.0 

KEL  -  (S-32.)*5/9  +  273.15 
H20HTVAPN  =  (2652.15+.6l6305*KEL-4.4788E- 
3*KEL**2)*7.7546 

DATA  F/l  .3088E8/,G/4.60E-l/.TCRIT/700.0/ 


TEGHTVAPN  -  (F*(1.-KEIVTCRIT)**G) 
♦4.3033E-4 

DATA  A/5.31 409/.B/-7.2001 7E-3/.C/  1 . 1 4153E-5/ 
HH20  =  V*(A*(KEL-273.15)  +B/2*(KEL**2- 
273.15**2) 

*  +C/3*(KEL**3-273.15**3))*7.7546  + 

V*H20TVAPN 
DATA  D/1.538E5/.E/5.87E2/ 

HTEG  =  U*(D*(KEL-273.15)  +  E/2*(KEL**2- 
273.15**2))*4.3033E-4 

*  +  U*TEGHTVAPN 
HVAP  -  HH20  +  HTEG 
RETURN 

END 


Figure  A-l.  Distillation  column. 
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Figure  A-2.  Equilibrium  tray  model. 


RTMENT  OF  THE  NAVY 

CM  Engineering  Laboratory 
ueneme.  CA  93043-5003 


l  Business 

/  for  Private  Use  $300 


Appendix  B 


ABSORBER  ANALYSIS 


The  absorber  column  is  modeled  using  the  equilibrium 
tray  approach  (see  Figure  B-l): 


Figure  B-l.  Absorber  column. 


0  -  Y„> 

3.  V  =  - V. - Mass  balance  on  air 

P  B 

(1  *  Y  ) 

4.  H,«f(Y1.V1,TUQ.B) 

H  =  f(Y  ,  V  ,  T  ) 

P  P  P  P 

5.  L  .  =  V  +  L_  -  V„  - Total  mass  balance 

6  X.,  =  (X.L,  +  Y  Vp  -  Y„  VB)/Lp. , - Water 

balance 

7.  hp, ,  -  (Lgh,  +  VH  -  VbHbVL  „ ,  - Energy 

balance 

«•  V.-KH^.L 
9.  Repeat  above  for  each  tray 


where:  V  *  total  molar  gas  flow  rate 

Y  =  mole  fraction  of  steam  in  gas 
H  -  inlet  vapor  enthalpy  for  gas 
L  -  total  molar  liquid  flow  rate 
X  -  mole  fraction  of  water  in  liquid 
h  -  enthalpy  of  water  in  liquid  at  inlet 

where:  p  -  any  plate 

B  -  the  bottom  of  the  column 
N  =  number  of  plates 

Known:  VB.  Y„.  H,.  L*.  X, 

Unknowns:  V  ,  Y  ,  H  ,  Lp. ,,  X  . hp, 1^,  N 

SOLUTION  PROCEDURE 

1.  Assume  T^  ^  \ 

2.  Y,  -  f(Xb)(equilibrium) 

Y  «f(X) 

o  'o' 


When  is  reached,  check  T^  against  known  value;  if 
necessary,  repeat  calculations  for  different  TLIQ  B 

PRCX3RAM  ABSORB 
REAL  LA  (100) 

REAL  LB 

DIMENSION  VA(IOO).  XA(100),  YA(100),  TO 00) 
DIMENSION  HAIR  (100),  HLIQ(IOO) 

KFLAG-0 
LB  -  25.2 
LA(1)=  LB 

PRINT,  “CHOOSE  ONE  OF  THE  FOLLOWING:" 

PRINT 

PRINT 

PRINT,  “1 -THERMAL  EQUILIBRIUM  (DEFAULT)" 
PRINT 

PRINT,  “2-ISOTHERMAL” 

PRINT 

INPUT,  TFLAG 
PRINT 

PRINT,  “ENTER  AIR  FLOW  RATE  IN  LB-MOLES/ 
HOUR:” 

INPUT.  VAIRIN 


B-1 


PRINT, 

PRINT,  “ENTER  EXIT  TEMPERATURE  OF 
LIQUID  IN  DEGREES  F” 

INPUT,  TBOT 
HLIQB  -  3645.6 

PRINT,  “ENTER  INLET  AIR  HUMIDITY  RATIO:” 
INPUT,  W 
YDUM  =  1.6092*W 
YAIRIN  =  YDUM/(  1  +  YDUM) 

PRINT,  “ENTER  INLET  AIR  TEMPERATURE  IN 
DEGREES  F:“ 

INPUT,  TAIRIN 
PRINT 

PRINT,  “ENTER  COLUMN  PRESSURE  IN 
ATMOSPHERES:" 

INPUT,  P 

YAIRIN  -  YAIRIN*P 

CALL  mixENTH( YAIRIN,  TAIRIN,  HAIRN) 

XT  -  .1445 
XEND  =  .020 
FINCON  -  .0001 
I-  1 

XB  *  .3182 
XA(1)  =  XB 
T(l)  =  TBOT 

50FORMAT(lX,TLATE’,’YA-,'XA’,’r,’V,T',’L’) 
PRINT  50 

100  CALL  EQUILA  (XA(I),t(i),YA(I)) 

YA(I)  -  YA(I)*P 

VA(I)  -  (VAIRIN  -  YAIRIN* VAIRINVO  -  YA(I)) 
CALL  mixENTH(YA(I),T(I),HAIR(I)) 

LA(I+1)  =  VA(I)  +  LB  -  VAIRIN 
XA(I+1)  =  (XB*LB  +  YA(I)*VA(I)  - 
YAIRIN*VAIRINyLA(I+l) 

HLIQ(I+1)  -  (LB*HLIQB  +  VA(I)*HAIR(I)  - 
V  AIRIN*HAIRINyLA(I+ 1 ) 

IF  (TFLAG.EQ.2)  THEN 
TO-1!  ) =  TBOT 
GOTO  110 
ELSE 

CONTINUE 
END  IF 

CALL  TLIQ(xa(i+l),HLIQ(I+l),T(I+l)) 

1 10  CONTINUE 

125  FORMAT  (lX.I3jF10.5,F9.5.2X,F6.2,F8.1.F8.1) 
PRINT  125, 1,  YA(I),  XA(I+1),  T(I),  VA(I),  LA(I+1) 
IF  (XA(I+1).LT.XT.AND.KFLAG.EQ.0)  THEN 
PLATES  -  I  -  1  +  (XA(I)-XTy(XA(I)-XA(I+l )) 
TTOP-  T(I)  +  (T(I+l)-T(I))*(XA(I)-XTy(XA(I)- 
XAO+l)) 

KFLAG-  1 
GOTO  130 
ELSE 

CONTINUE 


END  IF 

130  IF((XA(I)-XA(I+1).LT.0.0001)  GOTO  140 
IF  ((XA(I+l)-XEND).GT.FINCON)  THEN 
1  =  1  +  1 
GOTO  100 

ELSE 

CONTINUE 

ENDIF 

140  CONTINUE 
PRINT 

150  FORMAT  (IX, ’NUMBER  OF  THEORETICAL 
PLATES  =’E7.3) 

PRINT  150,  PLATES 

160  FORMAT  (IX, ’INLET  LIQUID  TEMPERATURE 
=’E7.3) 

PRINT  160,  TTOP 
170  END 
C 

subroutine  equila  (x.tsoln.y) 
real  xwatl(  1 3),xwat2(  1 1  ),slope(  1 1  ),yint(  1 3),bs(  1 1 ), 
cs(ll),  ds(ll), 

*  by(  1 3  ),cy(  1 3  ),dy(  13) 

real  x,w,tdew,tdewk,tsoln,pt,psat,m,b 
integer  n,  j 
character  flag*4 
c 

c  this  data  adapted  from  kohl  and  reisenfeld.  Figure  1 1.1 1, 
pg.  595 
c 

data  xwatl/0.0025,0.0042,0.0083, 0.0164,0.05,0.0776, 
0.10,0.2049, 

*0.3049,0.4808,0.6757,0.7813,0.8929/ 
dau  yint/-107.8.-105.3,-95.4,-88.2,-69.0,-61 .6.- 
56.0,-43.0,-34.3 

*  ,-24.5,-14.2,-9.0,-3.6/ 

dau  xwat2/0.0025 ,0.0042,0.0083 ,0.01 64,0.0402, 
0.02049,0.3049, 

*  0.4808,0.6757,0.7813,0.8929/ 

dau  slope/0. 4962,0.5 1 07,0.5 1 43 ,0.5393 .0.5654, 
0.7414,0.825, 

*  0.9091, 0.955, 0.9677.0.9882/ 
c 

pt=  101235. 

if  (flag.eq.'done’)  goto  2 
n  -  13 

call  spline  (n,  xwatl,  yint,  by,  cy,  dy) 

j-u 

call  spline  (j,  xwat2,  slope,  bs,  cs,  ds) 
flag  -  ’done’ 
c  alt  temp  in  deg  F 

2m  =  seval  (j,  x,  xwat2,  slope,  bs,  cs,  ds) 
b  -  seval  (n,  x,  xwatl ,  yint,  by,  cy,  dy) 
tdew  =  tsoln*m  +  b 
tdewk  -  (tdew  -  32.)*529.  +  273.15 


-2 


c  Maclaine-cross  1974  correlation  for  sat’n  pressure  of 
water  vapor 
c  all  pressures  in  Pa 

if  (tdewk  .It.  273.15)  then 

psat  -  exp<28.87  -  6133.8/tdewk) 
else 

psat  *  expf23.282-3780.82/tdewk-2258057 
tdewk**2) 

endif 

c  absolute  humidity  level  w  (mol  H20/mol  dry  air) 
w  «  0.622*psat/(pt-psat)*28.97/l8.02 
y  •  w/(l.  +  w) 

c  10  format  (f5.3,lx,f5.1,lx,f8.1,lx,f6.3,lx,f7.3,lx, 
f6.1,lx,f6.1, 

c  *  Ix,f7.1,lx,f7.5,lx,f5.3) 
c  write  (6,10)  x.tsoln,  ptjn.b.tdew, tdewk, psat, w,y 
return 
end 
C 

subroutine  tliq  (x.h.t) 

real  a,c,d,htsln,htsnmet,b,hmet,x,tk,t,aa,bb,cc 
c  A,B,C  are  coefficients  for  specific  heat  variation  of 
water  and  TEG 
c  resulting  units  are  J/kmol/K 

data  a/7J9e4/,c/lJ38e5/,d/5877jhtsln/1549.7/jef/ 
273.15/ 

hmet-b*1055/.454 
btsnmet-  hts  In  *  1 05 5/ . 45 4 

c  aa.bb.cc  are  all  coefficients  in  quadratic  formula 
aa-(1.0-x)*d/2.0 
bb-x*a  +  (1.0-x)*c 

cc--(x*a*ref  +  (1.0-x)*c*ref  +•  (1.0-x)*d/ 
2.0*ref**2  +  hmet 
*  +  x*htsnmet) 
c  units  of  TK  are  Kelvin 

tk«(-bb+sqrt(bb**2-4.*aa*cc))/27aa 
c  units  of  T  are  deg  F 
t«(tk-273.15)*975.+32.0 
c  10  format 

(fl0.1,lx,fl0.1,lx,l8.1>lx,fl0.1,lx,fl3.1,lx,f6.1,lx,f6.1) 
c  write  (6,10)  hmet,htsnmet,aa,bb,cc,ik,t 
return 
end 
C 

SUBROUTINE  MDCENTH  (y.T.hmix) 
c  Y  is  vapor  mole  fraction  of  water  in  air,  T  is  mixture 
temperature 

c  in  degrees  F;  HAIR,  HTVAPN,  HH20VAP,  and 
HMIX  all  in  BTU/lbmole. 
c  Assumption  regarding  temperature  of  mixture  is 
implemented  from  the 

c  main  program  by  passing  the  desired  temperature  to 
this  subroutine. 

REAL  y,TJcel,hair,htvapn,hh2ovap,hmix, 


a,b,c,d,e,f,g,h,i 
KEL  -(T-32.)*579.  +  273.15 
DATA  A/1.04195/, B/-2.63541E-4/,C/4.84486E=7/ 
HAIR  -  (l.-Y)*(A*(KEL-273.15)  +  B/2.* 
(KEL**2-273. 15**2) 

*  +C/3.*(KEL**3-273.15**3))*12.4667 
DATA  D/2652.15/ .E/.616305/E/-4.4788E-3/ 
HTVAPN  -  (A  +  B*KEL  +  C*KEL**2)*7.7546 
DATA  G/5.3 1 409/.H/-7.2001 7E-3/.I/1 . 1 41 53E-5/ 
HH20VAP-  Y*(G*KEL-273.15)  +  H/2.* 

(KEL**2-273. 15**2) 

*  +  I*(KEL**3-273.15**3))*7.7546  +Y*HTVAPN 
HMIX  -  HAIR  +  HH20VAP 

return 

end 

SUBROUTINE  SPLINE  (N,  X,  Y.  BY,  CY.  DY) 

INTERGER  N.NM1.IB.I 

REAL  X(N),Y (N),B Y(N),CY (N),DY(N) 

REALTY 

NM1-N-1 

IF  (N.LT.2)  RETURN 
IF  (N.LT.3)  GOTO  50 
DY(1)-X(2)-X(1) 

CY(2)-(Y(2)-Y(1VDY(1) 

DO  10I-2.NM1 
DY(I)-X(I+1)-X(I) 

BY(I)-2.*(DY(I-1)+DY(I)) 

CY(I+1)  -  (Y(I-t-l)-Y(I))/DY(I) 
CY(I)-CY(I+1)-CY(I) 

10  CONTINUE 
BYd)--DY(l) 

BY(N)  -  -DY(N-l) 

CY(I)  ■  0. 

CY(N)  -  0. 

IF  (N.EQ.3)  GOTO  15 

CY(1)  -  CY(3V(X(4)-X(2))  -  CY(2)/(X(3)-X(1)) 
CY(N)  -  CY(N-lV(X(N)-X(N-2))-CY(N-2)/(X(N- 
l)-X(N-3)) 

CY(1)-CY(1)*DY(1)**2/(X(4)-X(1)) 

CY(N)  -  -CY(N)*DY(N-l)**2/(X(N)-X(N-3)) 

15  DO  20  I-2JM 

TY-DY(I-1)/BY(I-1) 

BY(I)-  BY(I)-TY*DY(I-1) 

CY(I)  -  CY(I)-TY*CY(I-1 ) 

20  CONTINUE 

CY(N)  -  CY(NVBY(N) 

DO  30  IB-1J9M1 
I  -  N-IB 

CYd)  -  (CY(I)-DY(I)*CY(I+lVBYd) 

30  CONTINUE 

BY(N)  -  Y(N)-Y(NM1)VDY(NM1 )  + 

DY(NM1  )*(CY(NM  1  )+2.*CY(N)) 

DO  401-  1JSHV11 


B-3 


BY(I)  -  (Y(I+1)-Y(I)VDY(I)  - 
DY(I)*(CY(I+1  )+2.*CY(I)) 
DY(I)  -  (CY(I+1  )-CY (I)yDY(I) 
CY(I)  -  3.*CY(I) 

40  CONT1NIJE 
CY(N)  -  3.*CY(N) 

DY(N)  -  DY(N-l) 

RETURN 

50  BY(1)  =  (Y(2)-Y(l)y(X(2)-X(l)) 
CY(1) «  0. 

DY(1)  *=  0. 

BY(2)-  BY(1) 

CY(2)  -  0. 

DY(2)  «  0. 

RETURN 

END 


REAL  FUNCTION  SEVAL  (N,U,X,Y,B,C,D) 
INTERGERN 

REAL  U.  X(N),  Y(N),  B(N),  C(N),  D(N) 
INTERGER  I,  J,  K 
REAL  DX 
DATA  1/1/ 

IF  (I.GE.N)  I  =  1 
IF  (U.LT.X(I))  GOTO  10 
IF  (U.LE.XO+l)  GOTO  30 
10  I  -  1 
J-N+l 
20  K  -  (I+JV2 

IF  (U.LT.X(K))  J  =  K 
IF  (U.GE.X(K))  I  =  K 
IF  (J.GT.I+1 )  GOTO  20 
30  DX  =  U-X(I) 

SEVAL  -  Y(I)+DX*(B(I)+DX*(C(I)+DX*D(I))) 

RETURN 

END 


